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Abstract 

In the CP-violating Minimal Supersymmetric Standard Model (MSSM), for certain values 
of the CP-violating phases associated to the universal trilinear couplings {At,Af)) and the 
gluino mass (Mg), e.g., $cp = 60° or 90°, for Mh+ ^ 140 GeV and tan/? ~ 2 - 5, the 
lightest Higgs boson mass {Mhi) is ~ 50 GeV. This mass interval is still allowed from 
standard LEP Higgs searches because of a strongly suppressed HiZZ coupling. However, 
in the same region of parameter space in which these two conditions occur, the HiH^W^ 
coupling is enhanced because the two mentioned sets of couplings satisfy a sum rule. In 
this paper we probe such a light Higgs scenario at the Large Hadron Collider (LHC) by 
studying H^Hi associate production, leading to a 46 + £^ + I/!t signal. We show that the 
latter is readily accessible at the CERN hadron collider, upon the application of suitable 
selection cuts against the Standard Model (SM) backgrounds. Our parton level Monte 
Carlo (MC) analysis yields ~ 15 — 45 signal events, completely free of SM background, 
for C = 10 — 30 fb^^ of accumulated luminosity, after taking into account the overall 
efficiency for tagging four 6-jets. 
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The experimental observation of Higgs bosons and the determination of their properties is 
crucial for the understanding of Electro- Weak Symmetry Breaking (EWSB). Thus, the search 
for Higgs bosons is one the major goals of the present collider Tevatron (Run II) and future 
ones as well, such as the forthcoming LHC and the planned International Linear Collider (ILC). 
In the SM, the Higgs boson mass is not predicted. Over the last few decades, many efforts have 
been put into detecting such a particle, but to no avail. From direct searches at LEP, a lower 
bound of 114 GeV has been set on its mass |I1I2|- In the MSSM, with all real and CP-conserving 
parameters, the lower limit on the lightest Higgs boson is ~ 90 GeV [S] for any tan j3. However, 
this bound is significantly lowered in an MSSM scenario with radiatively induced Higgs sector 
CP- violation, as the latter in turn implies a suppressed HiZZ coupling j^. 

CP- violation in the Higgs sector is possible in multi- Higgs doublet models, such as a general 
2-Higgs Doublet Model (2HDM) or indeed the MSSM. In the latter, it has been shown that, 
assuming universality of the gaugino masses (Mj, z = 1, 2, 3) at some high energy scale, the CP- 
violating MSSM Higgs sector can be parametrised in terms of two independent phases: that of 
the Higgsino mass parameter (also called /i term), i.e., Arg(/i), and that of the soft trilinear 
Supersymmetry (SUSY) breaking parameters, i.e., ATg{Af), with f = t,b. The experimental 
upper bounds on the Electric Dipole Moments (EDMs) of electrons and neutrons E] as well 
as of mercury atoms 7] may pose severe constraints on these phases. However, these limits are 
highly model dependent. In particular, it has been shown that one could still have large CP- 
violating phases yielding relevant EDM values all satisfying current experimental bounds if any 
of these three possibilities is realised: (a) the sfermions of the first two generations are heavy, 
of the order of a few TeV 0; (6) cancellations between different EDM contributions [^j take 
place; (c) universality of the trilinear scalar couplings Aj is dismissed jTUllIII. However, only in 
the scenario with first and second generation sfermions much heavier than the third generation 
ones, the phase of /i can be large. Otherwise, it is strongly constrained, as Arg(/i) ^ 10~^. 

In the MSSM, non-zero phases of /i and/or the trilinear scalar couplings Af can induce 
CP-violation at one-loop level in the Higgs sector even in presence of a CP-conserving tree- 
level scalar potential, through the CP-violating interactions among Higgs bosons and heavy 
sfermions. This one-loop corrected Higgs potential then generates non-zero off-diagonal terms 
M |p in the 3x3 neutral Higgs boson mass-squared matrix Mfj, representing mixing between 
CP-even (or scalar, S) and CP-odd (or pseudoscalar, P) states These off-diagonal 

entries can approximately be written as follows 
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where "I'cp 
be 



Arg(y4f) = Arg(yu) and v = 246 GeV. The mass scale M|usy is typically defined to 
i.e., in terms of the two stop masses. After diagonalising the 3x3 symmetric 
Higgs mass-squared matrix Aijj by an orthogonal matrix O, the physical mass eigenstates 
Hi, H2 and (in ascending order of mass) are states of indefinite CP-parity. In this case, 
Mh± is the most appropriate mass parameter to describe the MSSM Higgs-sector (in place of 
Ma, often used in the CP-conserving case). 
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The CP-violating phases can cause the Higgs couphngs to fermions and gauge bosons to 
change significantly from their tree-level values, with dramatic consequences for MSSM Higgs 
phenomenology at present and future colliders JHl- Here, we start by presenting the form of 
the MSSM couplings HiVV, HiHjZ and HiH~W~^ , where V = W^, Z, in presence of explicit 
CP- violation'^: 

Ch^vv = gmw j29mvv[H^W+W-^ + -^H,Z^Z% (2) 
i=l ^^w 

^H.H.z = ^ J2 9mH,z{H,d,H,)Z^, (3) 
Ch.h-w+ = ^ j29mH-w+ {H.d^H-)W+^, (4) 



1=1 



where 



gH.vv = OiiCosP + 02iSmP, (5) 

gmHjZ = Osi{cosp02j - sinpOij) - {i ^ j), (6) 

9h,h-w+ = 02i cos (3 - Ou sin /3 + iO^i. (7) 

These couplings obey the following sum rules: 

T.9kvv = 1, (8) 



i=l 



9HiVV+ I 9h,H'W+ ? = 1, (9) 
9HkVv = eijkgn.HjZ- (10) 

Hence, from the above formulae, one can see that - if two of the gutzz are known - then the 
whole set of couplings between neutral Higgs bosons and gauge bosons is determined. It is also 
interesting to see from Eq. (9) that - in the presence of large CP violating effects with large 
scalar-pseudoscalar mixing - a suppressed HiVV coupling means an enhanced strength of the 
HiH~W^ vertex, as intimated already. This enhancement will indeed play a significant role in 
our analysis. Equally important is the correlation between the mass of the charged Higgs state, 
Mh±, and that of the pseudoscalar state. Ma- In fact, a suppressed HiVV coupling implies a 
light pseudoscalar state which in turn leads to a light charged Higgs, in particular Mh± < rrit 
(the top quark mass). 

From Eq. (1), it is clear that sizeable scalar-pseudoscalar mixing is possible for a large value 
of the CP-violating phase $cp and/or of | yU | and | ^4^ |. In this respect, the CP- violating bench- 
mark scenario defined as CPX in |13| provides a suitable choice of MSSM parameters which 
maximises CP-violating effects and can then be used to study the most striking phenomeno- 
logical manifestations of explicit CP-violation in the MSSM Higgs sector. Such a scenario is 



•^Further details can be found in Ref. 
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defined as followed: 



Mq = = = MsusY, (11) 
^ = 4MsusY, \At\ = lAI = 2MsusY, (12) 
Arg(A) = Arg(A;,). (13) 

Recently, the OPAL Collaboration ^Hj reported their results for Higgs boson searches in a 
CP-violating MSSM Higgs sector using the parameters defined in the CPX scenario and found 
that for the CP- violating phases $cp = 60° and 90° and certain values of Mh± and tan (3, the 
lowest mass limit on the neutral Higgs is very light - at times even vanishing completely - hence 
resulting in some low Higgs mass windows in the M/^±-tan/3 plane which are still allowed by 
LEP data. In a nutshell, the reason for the existence of such regions is the fact that in the 
CPX scenario the lightest Higgs boson is almost CP-odd with highly suppressed couplings to 
ZZ pairs. 

The experimental analysis was done by adopting both CPSuperH [20] and FeynHiggs 2.0 [21], 
as these are the two public codes available for the calculation of masses and mixing angles in 
the CP-violating MSSM Higgs sector. In reality, these two programs give somewhat different 
results - at least in the case of the CPX scenario - mainly due to different approximations used 
in their calculations. To give more conservative constraints, the experimentalists used the lower 
prediction of the two for the expected Higgs boson cross-sections. The constraints also depend 
sensitively on the mass of the top quark used in the calculation The results of jS], from a 
combined analysis of all LEP data, provide exclusion regions in the Mh^^ — tan f3 plane for the 
following values of the SUSY parameters: 



Arg(A) = Arg(A,) = ATg{M-,) = $cp, (14) 

MsusY = 0.5 TeV, Mg = 1 TeV, (15) 

M^ = My^ = 0.2 TeV, (16) 

$cp = 0°,30°,60°,90°, (17) 



where M^,M^^ and Mg represent the soft SUSY-breaking masses for bino, wino and gluino. 

By combining the results of Higgs searches from ALEPH, DELPHI, L3 and OPAL, the 
authors of Ref. [22] also provided exclusion regions in the Mh^ — tan (3 (as well as Mh+ — tan /?) 
plane for the same set of parameters. The exact shape of the exclusion regions may be somewhat 
different in their analyses, but they all show that for certain values of CP-violating phases 
LEP cannot rule out a light Higgs mass at low values of tan/5. This interesting situation 
roughly corresponds to tan/? ~ 3.5 — 5, Mh± ~ 125 — 140 GeV (yielding Mh^ ~ 50 GeV) and 
tan /3 ~ 2 - 3, Mh± ~ 105 - 130 GeV (yielding Mh, ^ 40 GeV), respectively. The authors of 
Ref. j22| further showed that in the same regions the Hitt coupling is suppressed too. Thus, 
these two particular areas of the parameter space can be probed neither at the Tevatron - 
where the WHi associated production mode is the most promising one - nor at the LHC - as 

''They were obtained for mt = 175 GeV, the value we adopt here. 
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the reduced Hitt coupling suppresses both the inclusive production mode gg — > Hi and the 
associated one tiHi. 

We have however found that, in the mentioned regions of the Mh+ — tan/3 plane, the de- 
cay — * HiW^ has a very large (~ 100%) Branching Ratio (BR), thanks to the discussed 
enhancement of the H'^HiW^ coupling and the mass hierarchy rrit > Mh± ^ Mh^- This 
feature motivated us to then study the possibility of probing such a light Higgs scenario in the 
CP-violating MSSM Higgs model through the process pp HiH^ — > Hi{HiW^) — > bbbblv at 
the LHC Recently, in Ref. j2Sl5 the authors have probed this light Higgs scenario through 
ti production at the LHC, where one of the top quarks decays into bbbW, via the chain 
t bH^, W^Hi and Hi — > bb, leading to a Ab + jj + + signal. In our analysis, the 

signal will consist of up to four 6-jets along with a hard lepton (electron or muon) and missing 
transverse energy J^t, according to the following decay pattern: 

PP^ H^ + Hi + X 




^ £u ^ bb 



In short, such a production and decay mechanism should allow to probe the possible existence 
of a light Hi state in the above mentioned two interesting windows of the CP-violating MSSM 
in the CPX scenario. The production mechanisms pp H^Hi and pp H'^H^ in the CP- 
violating MSSM and 2HDM, respectively, were discussed in Ref. where it was remarked 
upon the large if^ HiW^ BR, yet no phenomenological analyses were reported there. 

Our numerical results are obtained from a parton level MC analysis, wherein partons are 
treated as jets. As acceptance and selection criteria we required: 

1. I ?7 |< 2.5 for all jets and leptons, where rj denotes pseudo-rapidity; 

2. p^-j'^*^ > 15 GeV; 

3. pfp > 10 GeV (£ = e,/i); 

4. > 20 GeV; 

5. a minimum separation AR = \J {Acj))^ + {Arff = 0.4 between leptons and jets as well as 
each pair of jets; 

6. reconstruction of the two light Higgs bosons, by requiring four 6's in the event to be 
tagged as such^ and that at least one out of the three possible double pairings of 6-jets 

^With no jet- and/or lepton-charge determination, though. 
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satisfies the following mass constraint: 



where am = 0.12 Mhi- 

(In enforcing the latter constraint, we implicitly assume that trial values for Mh-^ are attempted 
and the selection has hit on the right one.) Notice that we impose Gaussian smearing on 
energies, with ^E/E = Q.Q/^/E for jets and AE / E = 0.12/Ve for leptons, to realistically 
emulate finite experimental resolution. All such cuts and smearing procedure have been applied 
to any process studied in this paper, alongside using CTEQ(5L) [SHI as Parton Distribution 
Functions (PDFs) taken at the scale = s (same for the scale of as, where relevant). 

In Figure 1 we show the variation of the signal cross-section (including the suppression due 
to a quadruple 6-tagging efficiency, = (1/2)^ = 1/16) with Mh+ [{a) and (6)] and Mh^ 
[(c) and (d)] for the CP- violating phase choices $cp = 60° [(a) and (c)] and 90° [(6) and (d)], 
respectively, for three values of tan f3. The choice of other MSSM parameters is defined through 
Eqs. (14)-(16). We have used the CPSuperH program |20J with rrit = 175 GeV to calculate the 
masses, couplings and decay rates of the relevant Higgs bosons and semi-analytical techniques 
to evaluate the hadro-production cross-section and decay rates. From Figure 1(a) and (b) one 
can see that the signal cross-section has a peculiar dependence upon Mh±. This may seem 
counter-intuitive, as the light Higgs mass increases with increasing Mu+. However, it should 
be noticed that, at lower Mh^ values, the &-jets emerging from the Higgs decays are rather soft 
and close to each other in phase space. As the light Higgs mass increases though, 6-jets become 
harder and also acquire much larger angular separations, hence a kinematics satisfying the cuts 
mentioned above more often, counter-balancing the decline in production rates due to larger 
values, ultimately leading to an overall relative rise in the cross-section. The final drop 
in the latter is due to phase space suppression for the H"^ — > HiW^ decay. In this scenario 
the largest cross-section, ~ 1.36 fb, can be obtained for $cp = 60°, tan/3 = 2 and Mh+ ~ 130 
GeV, which corresponds to Mh^ ~ 40 GeV. For the CP violating phase $cp = 90°, the largest 
cross-section, ~ 1.5 fb, can be obtained for tan/3 = 4 and Mh+ ~ 139 GeV, which corresponds 
to Mh^ ~ 50 GeV. From Figure 1(6) and (d), it is interesting to notice that the cross-section is 
almost insensitive to tan (3 at low values of Mh^ ■ 

The SM background cross-section arising from^ 

1. QCD production of gg —* bbjjii'; 

2. Electro- Weak (EW) qq' —* ZZW^ production, followed by the decay of each Z into bb 
pairs and by electron/muon decays of the VT-boson; 

3. top-quark production and decay via gg ^ tt ^ bbjjii'; 
^Hereafter, j labels both light {u,d,s,c) and heavy (6) quark jets. 
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is not shown, as it is negligible. In fact, after applying the same cuts 1.-6. to both signal and 
background processes and folding the cross-sections with the usual 6-tagging efficiency (ej, = 1/2 
per each &-jet) and the appropriate light- quark jet rejection factors (e.g., assuming Ru,d,s = 1/50 
and Rc = 1/25) 25j, we found that 

1. a{gg bbjjiu) < 2.2 x lO^^ fb; 

2. a{qq' ZZW^ bbjjiu) < 4.0 x lO'^ fb; 

3. a{gg ^ tt ^ bbjjiu) < 2.9 x 10"^ fb; 

where we have taken into account all relevant Z, BRs and combinatorial factors. 

In Ref. 1221 , the authors discussed qq' — > H2W^ as another possible probe of CP-violation 
in the MSSM Higgs sector, wherein iu and H2 — > HiHi bbbb. This mode can then 

lead to the same signature as the one we have been considering. Hence, one should in principle 
worry about its numerical relevance as well as possible interference effects between the two 
channels. In practise though, the qq' — *• H2W^ —>■ HiHiW^ — > bbbbiu production and decay 
rates are much smaller than those for qq' HiH"^ —>■ HiHiW^ —>■ bbbbiu. This is clearly 
confirmed by Figure 2. This is also the case for the interference. Therefore, hereafter, we will 
neglect considering the qq' —>■ H2W^ channel further. 

Note that signal events will be very striking due to the clustering of two pairs of bb invariant 
masses around M^^- This feature will of course not be present in the backgrounds, so it can 
be used to enhance the former and suppress the latter. One can attempt to reconstruct the 
light Higgs mass in the following ways. Out of the 46 final state, one can simply plot all 
possible (six) combinations of invariant masses (each with identical weight). This leads to 
the signal and background 'average' m^g distributions appearing in Figure 3(a). Alternatively, 
one can construct the three possible double pairings of bb invariant masses, then select the 
pairing giving the least difference between the two m^^ values, and the best reconstructed Hi 
mass is the corresponding mean value of that pair. This 'best reconstructed' distribution 
is shown in Figure 3(6). By comparing the spectra in Figures 3(a) and 3(6), it is clear that 
the second procedure is more efficient than the first one in increasing the signal-to-background 
rate. (Hence, we have exploited it by adopting the mass constraint described in Eq. (fT5|) .) For 
the two discussed mass spectra, two sample values of the light Higgs boson mass are assumed 
in the signal. Also note the normalisation to unity for all processes considered. Naturally, in 
producing these plots, we have refrained from applying cut 6. above. Figures 3(a) and 3(6) 
clearly highlight the Hi resonant peaks for the signal and the Z ones for the EW background. As 
for the other two background processes, the QCD one reveals a typical Jacobian shape, owning 
to the absence of heavy particles decaying hadronically, while the ti one displays resonance 
effects, when two light quarks are mistagged as 6-jets. In all cases, it is worth noting that the 
Gaussian smearing we have applied to the momenta somewhat affect the actual locations of the 
resonant peaks (where relevant). Under any circumstances, it is clear from both figures that 
the low mass signal resonances are always located far away from the bulk of all background 
events appearing at high mass. 
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The question now arises of whether it is also possible to reconstruct the mass in the 
signal. Because of the presence of a neutrino escaping detection in the — >• £z/ decay, the 
actual charged Higgs resonance is not kinematically accessible. Besides, only one of the two hb 
pairs selected via Eq. (fTHj) actually comes from a —>■ HiW^ bbW^ decay chain. In order 
to obviate such potential problems, we have proceeded as follows. Firstly, by running our MC 
for the signal, we have verified that the two (hereafter, 'primary') 6's emerging from the decay 
of the Hi boson produced in the hard scattering in association with the charged Higgs boson 
have normally a higher energy than the two produced in the mentioned decay chain (hereafter, 
'secondary'). This is evident from Figure 4, obtained after all cuts 1.-6. have been enforced. 
Hence, it makes sense to identify, between the two bb pairs selected via Eq. (fTHj). the one with 
least total energy as the one produced in the mentioned decay chain. Secondly, we have 
defined the transverse mass, Mt, constructed from, on the one hand, the visible transverse 
momentum of the system formed by the 'secondary' bb pair plus the lepton and, on the other 
hand, the missing energy, i.e.. 



Mt = ]J (p't + pI + Pt + ^t) - (pt + P^ + Pt + I^t) ' (19) 

where quantities in boldface refer to three- vectors. Such a variable is sensitive, as clearly evident 
from Figure 5, to the underlying charged Higgs boson mass and thus can be used within a MC 
simulation to fit the latter. 

In summary, we have proved the feasibility of testing at the LHC light Higgs boson windows 
in the so-called CPX scenario of the CP-violating MSSM, wherein a Hi signal might have been 
lost at LEP due to a strongly suppressed HiZZ coupling. Specifically, we have concentrated our 
attention upon the following two MSSM parameter space regions: (i) 3.5 < tan /? < 5, Mh^ ^ 50 
GeV and (ii) 2 < tan /? < 3, Mh^ ^ 40 GeV, assuming a common CP-violating phase $cp = 90° 
and 60°, respectively. We have found that, in the above mentioned parameter space areas, a 
light charged Higgs boson {Mh± < frit) can decay to W^Hi pairs with large a branching fraction 
so that, combined with a sizeable H^Hi associate production rate, the yield of the emerging 
signature ib + + J^t is sufficient to isolate a CP-violating signal. About ~ 15 — 45 signal 
events, completely free of SM background, for £ = 10 — 30 fb~^ of accumulated luminosity, 
after taking into account the efficiency for tagging all four 6-jets, can be isolated. Furthermore, 
we have also discussed the possibility of measuring the light Higgs boson mass as well as the 
charged Higgs boson one. We expect that the parton level study presented in this paper will 
encourage the CMS and ATLAS collaborations to carry out further investigations of the MSSM 
in presence of explicit CP-violation in the Higgs sector. 
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Figure 1: Variation of the signal cross-section with Mh± [(a), (&)] and M^i [(c), (rf)] for the two 
values of CP- violating phases $cp = 60° [(a), (c)] and $cp = 90° [(6), (d)]. The choices of tan/3 
for each CP-violating phase are shown in the figure. 
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Figure 2: A comparison between the production rates of the two processes qq' H2W^ 
HiHiW^ bbbUu and qq' HiH^ HiHiW^ bbbbiiy, for $cp = 60° (a) and $cp = 
{b), for the three choices of tan/3 shown in the figure. 
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Figure 3: The 'average' (a) and 'best reconstructed' (6) m^^ distributions (as described in the 
text) for the signal (with Mh^ = 32.38 and 54.32 GeV, corresponding to Mh± = 129.55 and 
136.97 GeV, respectively, for tan/5 = 4.5 and $cp = 90°). The corresponding SM model 
background distributions are also shown. 
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Figure 4: The energy distributions of the two 'primary' and two 'secondary' 6-jets in the signal 
(as defined in the text), as obtained for the representative parameter choices shown inside the 
figure. 
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Figure 5: The transverse mass distribution (as defined in the text) for the signal, for Mh± = 
129.55 GeV (left panel) and Mh± = 135 GeV (right panel), tan/3 = 4.5 and $cp = 90°. 
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